TrES-5: A Massive Jupiter-sized Planet Transiting A Cool G-dwarf 

Georgi Mandushev 1 , Samuel N. Quinn 2 , Lars A. Buchhave 3 , Edward W. Dunham 1 , Markus 
Rabus 4 ' 5 , Brian Oetiker 6 , David W. Latham 2 , David Charbonneau 2 , Timothy M. Brown 7 , Juan 

A. Belmonte 5,8 , Francis T. O'Donovan 9 

ABSTRACT 

We report the discovery of TrES-5, a massive hot Jupiter that transits the star 
GSC 03949-00967 every 1.48 days. From spectroscopy of the star we estimate a stellar 
effective temperature of T e g = 5171 ± 36 K, and from high-precision B, R and I pho- 
tometry of the transit we constrain the ratio of the semi-major axis a and the stellar 
radius R± to be a/R* = 6.07±0.14. We compare these values to model stellar isochrones 
to obtain a stellar mass of M* = 0.893 ± 0.024 M@. Based on this estimate and the 
photometric time series, we constrain the stellar radius to be R+ = 0.866 ± 0.013 Rq 
and the planet radius to be R p = 1.209 ± 0.021 Rj. We model our radial- velocity data 
assuming a circular orbit and find a planetary mass of 1.778 ± 0.063 Mj. Our radial- 
velocity observations rule out line-bisector variations that would indicate a specious 
detection resulting from a blend of an eclipsing binary system. TrES-5 orbits one of the 
faintest stars with transiting planets found to date from the ground and demonstrates 
that precise photometry and followup spectroscopy are possible, albeit challenging, even 
for such faint stars. 

Subject headings: planetary systems — techniques: photometric — techniques: radial 
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Introduction 



Long before the initial di s coveries of extrasolar p lanets around Sun-like stars by the Doppler 



technique (jLatham et al. 



1989 



Mayor Queloall995l ). it was recognized tha t planets aro und other 



stars could also be detected photometrically if they transit their host stars (|Struvelll952l ). Transit 
observations can provide at least two crucial pieces of information: the planet's size and orbital 
inclination. When combined with radial velocity data they can be used to derive the planet's 
mass, surface gravity and density, and infer its internal structure and composition. Many other 
studies also become possible for transiting planets around bright host stars, including analysis of 
the planet's atmosphere, temperatu re distribution, the presen ce of other bodies in the system, the 



spin-orbit alignment, etc. (see, e.g.. ICharbonneau et al 



20071) 



The ever increasing number of discovered transiting planets, both from ground -based sur 



veys and from the Kepler and CoRo T missions (jBorucki et al 



2010 



Baglin et al.ll2009l ). shows the 



tremendous diversity of extrasolar planets and planetary systems. Particularly interesting are the 
relations between planetary mass, radius and orbital period, as well as the distribution of planetary 
masses, as together they hold the key to understanding planetary formation and evolution. The 
current summary of transiting exoplanets at http : / / exoplanet . eu indicates that 70% of the tran- 
siting planets with known masses have masses M p < 1.5 Mj. This apparent drop in the number 
of massive planets (noted also by ISouthworth et al.ll2009l ) cannot be explained by major selection 
effects. Massive planets are also large, and thus easier to discover, unlike lower -mass planets, whos e 
radius increases monotonically with mass for non-irradiated planets with cores (jFortnev et al.ll2007l ) . 



We report here the discovery of the Jupiter-sized massive transiting planet TrES-5. The planet 
orbits one of the faintest star yet (V = 13.72) with a transiting planet found from ground-based 
wide-field surveys, and is one of the very few planets found around stars fainter than V = 13.5, 
including the space-based Kepler and CoRoT missions. TrES-5 increases the relatively small number 
of well-characterized massive transiting plane ts. Such planets play an important role in the study 
of planet format ion, migration and evolu tion (jChatteriee et al.ll2008l ; iBaraffe et al.ll2010l ) and it has 
been suggested (jSouthworth et al.ll2009l ) that the massive extrasolar planets form a different planet 
population than the lower-mass ones. 



2. Photometry and Spectroscopy 



We monitored a 5? 8 x 5? 8 field in Cygnus with the Lowell Observatory Planet Search Survey 
Telescope (PSST, iDunham et al.ll2004l ) and the STARE telescope on the Canary Islands in Spain, 
between UT 2007 July 15 and UT 2007 October 7. Both telescopes are part of TrES (the other 
two telescopes — Sleuth, at Palomar Observatory, and WATTS, in Texas, were not operational 
at that time). The same field was observed again in 2009 and 2010 by P SST alone. All image s 
were reduced, and photometry and transit search carried out as described in lDunham et al.l ((2004). 
Altogether PSST observed 13 full and partial transit-like events of the star GSC 03949-00967. The 
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object was at the very edge of the PSST field and because of the slightly different field size and 
scale the star was not observed by the STARE telescope. The depth and duration of the events 
were consistent with the transit of a Jupiter-sized planet across a K dwarf, and we undertook a 
program of followup observations to confirm the planetary nature of the object and measure its 
properties. The discovery light curve of TrES-5 is shown on Fig. [1] 

We carried out high-precision in-transit B, R, and /-band photometry of TrES-5 using Lowell 
Observatory's 0.8m, 1.1m Hall and 1.8m Perkins telescopes. We observed TrES-5 on UT 2009 
November 17 with the 0.8m telescope with NASACam in Cousins R (47 exposures), on UT 2010 
June 5 and 2010 November 30 with the 1.1m Hall telescope with the NASA42 camera in Cousins 
R and I (87 and 108 exposures, respectively), and on UT 2010 September 5 with the 1.8m Perkins 
telescope with PRISM in Johnson B and Cousins I (64 and 57 exposures, respectively). For all 
data sets, we derived differential fluxes relative to a large ensemble of local comparison stars. The 
high-precision Lowell Observatory photometry is shown in Fig. [2j 

In order to characterize the host star and provide an external check on the stellar parameters 
derived from the TRES spectra, we obtained off-transit BV(RI)c photometry of the TrES-5 system 
on UT 2010 October 10 and 11 with the 1.05-m Hall telescope at Lowell Observatory in combination 
with a 4K x 4K e2v CCD231 CCD detector. We calibrated the photometry using 6 standard fields 



from ILandoltl (|1992l ) . The results are listed in Table [T] together with other relevant data for the 



host star of TrES-5. 



We observed the candidate with the Tillinghast Reflector Echelle Spectrograph (TRES; lFuresz 



2008), mounted on the 60- inch Tillinghast Reflector at Fred L. Whipple Observatory on Mt. 
Hopkins in Arizona from 2010 September to 2011 April. We obtained 8 spectra, each consist- 
ing of 51 echelle orders and spanning the wavelength range 3850-9100 A, with a resolving power 
of A/AA « 44,000. We obtained absolute radial velocities by cross-correlation against a syn- 
thetic template chosen from a large librar y of spectra based on Kurucz model atmospheres (see 



Nordstrom et al.lll994l : lLatham et al.ll2002l ). calculated by John Laird and based on a linelist com- 
piled by Jon Morse. These velocities have a typical precision of 0.2 kms" 1 . However, the synthetic 
spectra only span a small portion of the total wavelength range, centered on the Mg b triplet 
at 5187 A. We derived precise relative velocities by taking advantage of the full spectrum, cross- 
correlating all orders against those of an ob served template. For a description of the reduction and 



cross-correlation, see lBuchhave et al.l (|2010l ). in which the details of the same reduction package are 
described, as applied to the Fiber-fed Echelle Spectrograph (FIES). The absolute gamma velocity, 
calculated using the weighted mean offset of the relative velocities from the corresponding absolute 
velocities, was determined to be 7 = —13.403 ± 0.100 kms -1 . The absolute gamma velocity is onto 
the system in which HD 182488 is defined to have an absolute velocity of —21.508 kms -1 . 

To derive the stellar atmospheric parameters, the observed spectra were cross-correlated against 
a grid of synthetic spectra drawn from the library of spectra described above. The synthetic spectra 
cover a window of 300 A centered near the gravity-sensitive Mg b triplet and have a spacing of 
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250 K in effective temperature, 0.5 dex in gravity, 0.5 dex in metallicity and 1 kms -1 in rotational 
velocity. The best matched template to the observed spectrum represents the best matched stellar 
parameters on the library grid. A new set of tools is then used to derive more precise stellar pa- 
rameters from the normalized cross correlation peaks. A description of the tools will be published 
in Buchhave et al. (2011, in preparation). 

We used the multicolor Johnson- Cousins photometry to form a vari ety of color indices, and u sed 
the color-temperature calibrations bv lRamfrez &: Melendea (|2005l ) and lCasagrande et al.l (|2006l ) to 
estimate the star's effective temperature T e ff. The average result, T e g = 4943=L36 K, is 228 K cooler 
than the spectroscopic value of T e g- = 5171 ±36 K (see below). This result is not surprising as some 
reddening is expected at the object's Galactic latitude of 13°. Prom the match to stellar evolution 
models described below, we estimate a reddening of Eb-v ~ 0.07, which is in good agreement with 
the difference of 228 K between the spectroscopic and photometric effective temperatures. 



3. Light Curve Analysis 



We analyzed the five B, R and /-band transit observations using the analytical expressions in 
Mandel h Agoll (|2002l ) to compute the model flux. We assumed a circular orbit and a qu adratic 
stellar limb-darkening law, fixing the coefficients at the color-dependent values tabulated in IClaret 
(|2004l ) for the spectroscopically-estimated T e g, log <?, and [Fe/H]. Initially we adoped as free param- 
eters the orbital period P and epoch To, the ratio of the planet radius to the stellar radius -Rp/i?*, 
the square of the impact paramater b = (a/R+) cosi, where a is the semi-major axis of the planet's 
orbit and i is the orbital inclination, the parameter C/-^*> which for zero eccentricity is related to 
a/R± via (/R* = i a / '-R*)(27r y 'P)(l - 6 2 )" 1 / 2 , and the zero point of the out -of-transit flux F n . This 
choice of parameters minimizes the correlation between b and a/P* (see iBakos et al.l 120071 ). We 
performed an initial fit using the high-precision followup transit photometry, after which we fixed 
all parameters except P and To to their values from the initial fit, and repeated the fit using both 
the PSST wide-field photometry and the followup photometry. This approach allowed us to utilize 
the entire baseline of over three years of observations and obtain more precise values for P and To 
than if only the followup photometry were used. We then fixed the values of To and P (stated in 
Table [3]) in the subsequent analysis. 



W e performed a joint fit of all five transit light curves using the MPFIT package (jMarkwardt 
2008), and found the values of P p /P*, b 2 , (/R* and Fq that minimized the x 2 - This fit is shown 



with solid lines in Fig. (2) In order to estimate the errors of th e fitted parameters, we conducted a 
Markov Chain Monte Carlo (MCMC) analysis (see lFordll2005l . and references therein). Altogether 
10 independent Markov chains with 5 x 10 5 points per chain were created, each chain starting 
from a random perturbation of the best-fit parameters. We discarded the first 20% of the points 
to minimize the impact of the initial conditions. The resulting histograms of 4 x 10 6 values were 
used to obtain the mode and the 68.3% lower and upper confidence limits p\ Q and p up for each 
parameter. We then adopted the mode as the final, optimal value of the parameter, and the value 
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°f (Pup — Pio)/2 as its l-o" error (encompassing 68.3% of the parameter values around the mode). 
The use of a single number to characterize the width of each parameter's distribution is justified 
by the nearly Gaussian, symmetric shape of all histograms. 



4. Properties of TrES-5 And Its Star 



The mass and radius of the host star of TrES-5, required for establishing the planet properties, 
were determined on the basis of the spectroscopic T c g and [Fe/H], and the value of a/R* derived 
from the light curve fit described above. The quantity a/R+ is closely related to the stellar density, 
and is determined in thi s case with higher r elative precision than logg. It is therefore a better 
proxy for luminosity (see ISozzetti et al.ll2007l ). 



m 



We used an iterat ive procedure to determine the star's parameters, similar to the on e descr ibed 



Torres et al.l (|2008l ). Using the stellar evolution models from the series bv lYi et al.l (|200ll ). we 



computed isochrones over a range in stellar age, from 0.1 Gy to 12 Gy, over the metallicity range 
allowed by the spectroscopy. We compared the measured T e g and a/R* from the light curve fit to 
the computed values of a/R* and T e s along each isochrone and recorded the points on the isochrone 
which matched the measured quantities within their errors. We then calculated the weighted mean 
of all matches, with weights inversely proportional to the exponent of the distance (in \ 2 sense) 
between the observed and model values of T e g, a/R+ and [Fe/H] (see iTorres et al.l 120081 ) . 



This procedure yields a value of log g which is better constrained than the spectroscopic es- 
timate. With the value of log g fixed, revised estimates of T e g and [Fe/H] were derived from the 
spectra, and the light curve fit repeated. The new values of T e g and [Fe/H] (from spectroscopy), 
and a/R* (from the fit) were again compared to the stellar models to obtain the final values of 
Af* = 0.893 ± 0.024 M , R* = 0.866 ± 0.013 R & , logp = 4.513 ± 0.013, and age of 7.38 ± 1.87 Gy. 

We fit a Keplerian orbit to these data assuming zero eccentricity as a good first approximation, 
as expected from theoretical arguments for a period as short as 1.48 days. The period and epoch 
were held fixed. The RMS of this fit is 24.2 ms" 1 , which is similar to the internal errors of the 
velocities. The parameters of this orbital solution are listed in Table El The orbit is displayed in 
Fig. [3] (top panel) along with the observations, and the residuals are shown in the middle panel. 

We investigated the possibility that the radial velocities we measured are the res ult of distor 



tions in the line profiles due to conta mination from an un resolved eclipsing binary (j Santos et al 



2002 



Torres et alJl200a ) or star spots (IToner Graylll988l ). instead of being due to true Doppler 



motion in response to a planetary companion. We cross-correlated each TRES spectrum against 
a synthetic template matching the properties of the star, and averaged the correlation functions 
over the orders. From this representation of the average spectral line profile we computed the mean 
bisectors, and as a measure of the line asymmetry we calculated the "bisector spans" as the velocity 
difference between point s selected near the top and bottom of the mean bisectors (jToner fc Gray 



1988 



Torres et al.l 120051 ) . If the velocities were the result of a blend with an eclipsing binary 
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(jMandushev et al.ll2005l ) or star spots ([Queloz et al.ll200ll ). we would expect the line bisectors to 
vary in phase with the photometric period with an amplitude similar to that of the velocities. In- 
stead, we detect no variation in excess of the measurement uncertainties (see Fig. O bottom panel), 
and we conclude that the velocity variations are real and that the star is orbited by a Jovian planet. 



Discussion 



TrES-5 orbits one of the faintest host stars with transiting planets found to date. The discovery 
demonstrates that precise photometry and spectroscopy are possible from the ground even for such 
faint stars, although some followup studies will be more challenging. With a period of only 1.5 days, 
TrES-5 is a classical "hot Jupiter", although the relatively cool host star (T e g = 5171 K) places 
it near the upper (hotter) end o f the transition zone between the "pM class" and "pL class" 
planets proposed in iFortnev et al.l ()2008l ). In such planets the importance of gaseous TiO and VO 
absorption is diminished, temperature inversion may not be present, and the flux from the host 
star is more evenly distributed around the planet and inside its interior. Because of its location 
near the pM/pL boundary, TrES-5 could be a suitable object for testing the models of planetary 
atmospheres. 

At M p = 1.8 Mj TrES-5 is above the upper quartile of the planetary mass distribution. The 
mass and radius o f TrES-5 are in good agreement with the th eoretical mass-radius relations for 
irradiated planets (jFortnev et al.l 120071 : lHansen fe Barman! 120071 ). For an age of 7.5 Gyr, planetary 



mass of M p = 1.8 Mj and equilibrium temperature of T eq ~ 1500 K, the predicted radius of the 
planet is R p ~ 1.2 Rj, which is in good agreement with the observed value. Thus, TrES-5 does 
not appear to have the anom alously large radius of many close-in giant exoplanets. A comparison 
with the planetary models of IFortnev et al.l (|2007l ) for the mass, age and semimajor axis of TrES-5 
indicates that the planet is best approximated by a model with no heavy-elements core. 

It has been noted t hat many massive planets have non-zero eccentricity even at short periods 
(jSouthworth et al.ll2009l ). Using the current catalog of transiting planets at http : / / exoplanet . eu, 
we find that over a third of the short-period, massive planets (P < 3 days and M p > 1.7 Mj) have 
non-zero eccentricities. In our radial velocity fit we have assumed e = based on the short period 
of Tr ES-5 and theoretical arguments about the time scales of tidal interactions (IZahn Bouchet 
19891 ). Our radial velocities clearly reject large eccentricities. From MCMC analysis of the radial 
velocity data, we find that the prefered orbital solution has e = 0.0251q q25- We therefore adopt 
e = for the orbit of TrES-5, but we cannot rule out a small eccentricity (e < 0.04). A future 
detection of the secondary eclipse will help constrain the orbital eccentricity much better, as well 
as test for the presence of temperature inversion in the atmosphere of TrES-5. 
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Fig. 1. — The discovery light curve of TrES-5. The plot shows the relative flux of the TrE5-5 
system as a function of the orbital phase, adopting the ephemeris in Table O 
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Fig. 2. — High-precision followup B, V and ii-band photometry of TrES-5. The plot shows the 
relative flux of the TrE5-5 system as a function of time relative to the center of transit, adopting 
the ephemeris in Table [3j Each light curve is labeled with the telescope and date of observation. 
The residuals from the simultaneous fits (overplotted with solid lines) are shown below each light 



curve. 
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Fig. 3. — Top: Radial velocity observations of TrES-5 obtained with TRES, shown relative to the 
center of mass and adopting the ephemeris in Table The best-fit orbit (solid line) is overplotted. 
Middle: Residuals from the best-fit model to the radial velocities. Bottom: Bisector spans shifted 
to a median of zero, for each of the TRES exposures. 
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Table 1. TrES-5 Host Star 



Parameter Units Value Source 



RA 


J2000.0 


20 h 20 m 53 ! ?24 


1 


Decl. 


J2000.0 


+59°26'55'. / 6 


1 


GSC 




03949-00967 




r 1 

[Ma, m 


_i 

mas yr 


[+10.4, +30.5] 


1 


V 




13.718 + 0.005 


2 


B-V 




0.927 + 0.008 


2 


V-R c 




0.526 ± 0.007 


2 


V-I C 




0.994 ± 0.007 


2 


J 




12.111 + 0.027 


3 


J-H 




0.438 ± 0.035 


3 


J-K s 




0.547 + 0.033 


3 




M & 


0.893 ± 0.024 


2 


R* 


Rq 


0.866 ± 0.013 


2 


T c s 


K 


5171 ± 36 


2 


[Fe/H] 




+0.20 + 0.08 


2 


logs' 




4.513 + 0.013 


2 


vsini 


kms -1 


3.80 ± 0.36 


2 


My 




+5.756 + 0.041 


2 


Distance 


pc 


360 + 11 


2 


Age 


Gy 


7.38 ± 1.87 


2 



References. — (1) UCAC 3 (IZacharias et al.ll20im : 
(2) this paper; (3) 2MASS (|skrutskie et alJbood ) 
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Table 2. Radial velocity measurements of TrES-5 



HJD RV (ms- 1 ) 



2455456.847578 


-542.8 


± 


27.9 


2455468.707361 


-472.6 


± 


23.3 


2455647.005138 


-308.3 


± 


22.8 


2455649.002656 


127.2 


± 


26.9 


2455651.985839 


98.3 


± 


29.0 


2455659.987762 


-549.8 


± 


29.9 


2455667.970142 


0.0 


± 


14.8 


2455673.945954 


45.0 


± 


14.8 



Table 3. TrES-5 Planet Parameters 



Parameter 


Units Value 


P 


days 1.4822446 ± 0.0000007 


To 


HJD 2 455 443.25153 ± 0.0001069 


a 


AU 0.02446 ± 0.00068 


i 


deg 84.529 ± 0.005 


a/R* 


6.074 ±0.143 


b = acosi/R* 


0.579 ± 0.026 


K 


ms" 1 339.8 ±10.4 


M p 


Mj 1.778 ±0.063 


Rp 


,Rj 1.209 ±0.021 


P 


gem" 3 1.25 ±0.08 


T 


K 1484 ± 41 


Rp/ R* 


0.1436 ±0.0012 



Note. — Mj = 1.899 x 10 27 kg (Jupiter's mass); R.i = 
7.1492 x 10 7 m (Jupiter's equatorial radius) (|Coxll200ol ) 



